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Nonr-4104  (00). 

The  program  was  sponsored  by  the  Advanced  Research  and  Project  Agency  through  the  Power 
Branch  of  the  Office  of  Naval  Research  under  the  direction  of  Dr.  J.  Huth  of  ARPA  and  Mr. 

J.  A.  Satkowski  of  ONR. 
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I.  INTRODUCTION 


This  second  Quarterly  Technical  Summary  Report  describes  the  progress  made  in  the  period 
1  August  through  October  31.  During  this  period,  the  construction  of  the  Allison  diagnostic 
MPD  device  was  completed  and  the  first  checkout  run  was  made.  This  report,  therefore,  deals 
with  the  final  description  of  the  loop.  Most  emphasis  is  laid  upon  the  description  of  the  diagnostic 
methods,  which  were  prepared  for  the  use  on  the  loop.  Therefore,  actual  test  data  with  the 
unseeded  working  fluid  will  be  included  in  the  next  quarterly  report. 
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II.  RESUME  OF  PROGRESS 


CONSTRUCTION  OF  CLOSED  LOOP  DIAGNOSTIC  MPD  DEVICE 

The  final  assembly  of  the  diagnostic  MPD  system  has  been  completed.  For  the  first  checkout, 
the  tantalum  MPD  test  section  was  replaced  by  a  stainless  steel  insert.  The  tantalum  MPD 
test  section  itself  has  been  completed  and  is  ready  for  incorporation  into  the  system. 


The  system  was  operated  for  a  total  of  one  hour  and  then  taken  apart  for  inspection  of  the  com¬ 
ponents.  The  run  with  the  highest  temperature  was  made  at  1000'’K.  The  impurities  measured 
were  oxygen  with  1.  75  ppm  and  water  vapor  with  48  ppm. 

A  sample  piece  of  tantalum  was  weighed  and  then  installed  in  the  stainless  steel  insert  section 
for  the  aforementioned  test.  It  did  not  show  any  wear  after  the  test;  therefore,  more  running 
time  must  be  accumulated  before  the  weight  of  the  lost  mass  can  be  determined. 

A  gas  discharge  tube  was  manufactured  which  can  be  attached  to  the  insert  test  section  for 
spectroscopic  determination  of  the  impurities.  This  investigation  is  scheduled  for  the  next  run 

Presently,  it  is  planned  to  accumulate  more  running  time  on  the  system  to  get  more  informa¬ 
tion  on  the  reliability  of  the  single  components.  *  Following  this,  the  tantalum  test  section  will 
be  incorporated  into  the  system  and  the  measurements  with  the  unseeded  working  fluid  will  be 
completed  in  the  next  reporting  period. 

DEVELOPMENT  AND  ADAPTION  OF  DIAGNOSTIC  METHODS 

Further  development  of  diagnostic  methods  was  accomplished.  A  part  of  this  work  was  the 
subject  of  an  article,  "Calculation  and  Measurement  of  the  Radial  Distribution  of  Line  Broad¬ 
ening  in  a  Cesium  Plasma,  "  by  D.  L.  Tipton  and  R.  T.  Schneider,  Applied  Spectroscopy  17. 
147,  1963.  This  work  is  an  application  of  Abel's  integral  equation  on  the  broadening  of  a 
spectrum  line  to  determine  the  electron  density  profile  of  a  cylindrical  plasma.  This  method 
will  be  applied  to  find  changes  in  the  electron  density  in  the  MPD  device. 


=!Tn  the  meantime,  eight  hours  of  running  time  have  been  accumulated  at  1350°K  maximum  tem¬ 
perature. 


Allison 


Another  part  of  the  work  accomplished  during  the  reporting  period  is  described  in  Spectro¬ 
scopic  Technique  for  Temperature  or  Pressure  Measurement  in  Plasmas,  by  R,  T,  Schneider, 
W,  Woerner,  and  H,  E,  Wilhelm,  This  paper  included  as  a  part  of  Section  IV  of  this  report, 
describes  a  new  method  for  spectroscopic  temperature  measurement.  The  advantages  of  the 
new  method  are:  (1)  the  transition  probabilities  need  not  be  known,  and  (2)  the  photographic 
plate  need  not  be  calibrated, 

THEORETICAL  INVESTIGATIONS 

Ionization  Equilibrium  in  Magnetic  Field 


The  partition  functions  for  the  components  of  a  partially  ionized  plasma  energized  with  a  mag¬ 
netic  field  have  been  calculated  on  a  statistical  basis.  These  allow  a  complete  description  of 
the  magnetothermodynamic  properties  of  the  plasma.  As  an  application,  the  ionization  equili¬ 
brium  in  a  strong  magnetic  field  will  be  calculated. 

Nonequilibrium  Ionization  in  Magnetic  Field 


Statistical  investigations  of  the  properties  of  a  reacting  plasma  by  means  of  the  kinetic  equations 
have  been  started.  As  an  intermediate  application,  the  electrical  and  thermal  conductivity  of 
seeded  plasma  will  be  calculated.  Following  this,  a  kinetic  nonequilibrium  ionization  theory 
will  be  developed 

Results  of  these  investigations  will  be  included  in  the  third  quarterly  report. 
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III.  THE  ALLISON  DIAGNOSTIC  CLOSED  LOOP  MPD  DEVICE 


FINAL  ASSEMBLY  AND  LEAK  TESTING 

The  design  of  the  Allison  MPD  diagnostic  loop  was  described  in  the  first  Quarterly  Report, 

EDR  3511.  Therefore,  in  this  report  only  the  final  assembly  is  described  and  a  few  minor 
changes  in  the  design  are  pointed  out.  Figure  1  is  a  schematic  of  the  loop.  The  changes  in  this 
schematic  as  compared  with  the  schematic  presented  in  EDR  3511  are: 

•  An  additional  line  was  installed  to  bypass  the  helium  compressor.  The  reason:  to  take  the 
compressor  (a  possible  major  leak  source  during  the  mass-spectrometer  leak  testing)  out 
of  the  system.  Once  the  system  was  made  leakproof,  the  compressor  was  put  back  in  the 
system;  any  new  leaks  now  could  be  located  more  easily. 

•  A  new  vacuum  line  was  installed  for  evacuating  the  bearing  housings  separately.  The 
narrow  clearance  precluded  a  good  evacuation  through  the  main  vacuum  line.  Also,  the 
purified  helium  refill  will  be  made  with  the  same  line.  Another  sampling  valve  was  in¬ 
stalled  for  continuous  sampling  in  the  vicinity  of  the  orifice  flow  meter. 

•  A  bypass  line  from  the  MPD  test  section  back  to  the  system  was  installed  to  handle  the 
helium  cover  gas  which  is  used  to  cool  the  test  section. 

Figure  2  depicts  the  system  during  assembly;  this  is  a  portion  of  that  part  of  the  system  which 
is  located  inside  the  building  (also,  see  Figure  1).  The  following  components  are  shown  in 
Figure  2:  charcoal  filter,  NaK  bubblers,  magnet,  cesium  tank,  helium  bottles. 

The  part  of  the  system  outside  the  building  is  shown  in  Figure  3.  Included  in  Figure  3  are  the 
following  components;  helium  compressor,  filter,  cesium  tank,  helium  bottles. 

Figure  4  is  a  special  view  of  the  coolers.  The  coolers  are  sloped  so  that  the  cesium  can  flow 
back  to  the  helium  cesium  separator  located  at  the  extreme  left  in  Figure  4. 

Figure  5  shows  the  control  panel.  The  location  of  the  controls  which  are  panel  mounted  is 
shown  on  the  schematic  in  Figure  1. 

The  heater  and  the  magnet  are  depicted  in  Figure  6.  The  hole  in  the  core  of  the  magnet,  which 
gives  access  for  optical  observation  of  the  test  section,  is  visible  in  the  center  of  Figure  6. 
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After  assembly,  the  system  was  leak  tested  using  a  mass  spectrometer  leak  detector.  Only 
a  few  small  leaks  occurred  in  the  system  itself.  None  of  the  welds  had  to  be  redone.  Approx¬ 
imately  ten  medium  leaks  had  to  be  repaired  in  the  helium  compressor.  After  this  was  done, 
the  system,  with  the  helium  compressor  and  a  stainless  steel  insert  test  section  in  place,  could 
be  pumped  down  to  1  X  10"^  Torr. 

FIRST  CHECKOUT  RUNS 

The  main  purposes  of  the  first  checkout  runs  were  to  test  the  heater  and  to  get  preliminary 
information  on  the  impurity  level  in  the  system. 

Figure  7  shows  the  heater  coil  before  the  first  test.  Figure  8  shows  the  same  coil  after  one 
hour  running  time  with  a  maximum  gas  temperature  of  1000‘’K.  The  tungsten  wire  has  been 
polished  by  the  high  velocity  helium  gas  at  those  areas  where  a  high  temperature  was  encoun¬ 
tered.  See  Figure  8.  A  part  of  the  coil  looks  bright  metallic,  while  before  the  run  the  surface 
of  the  coil  was  dark  and  rough  (see  Figure  7).  Traces  of  tungsten  were  evaporated  from  the 
coil  and  condensed  on  the  insulator.  This  also  can  be  seen  when  comparing  the  white  insulator, 
which  supports  the  tungsten  coil  in  Figure  8,  with  that  in  Figure  7. 

For  this  preliminary  run  a  boron  nitride  insulator  was  used.  The  boron  nitride  insulator 
worked  out  better  than  expected.  However,  for  higher  gas  temperatures  thoria  will  be  used  to 
avoid  |x)ssible  damage  of  the  tantalum  test  section  by  nitrogen  and  of  the  stainless  steel  coolers 
by  boron.  ThOria  insulators  are  on  order  and  will  be  available  for  the  actual  MPD  runs.  The 
boron  nitride  insulator  shows  some  damages.  Small  pieces  broke  off,  probably  due  to  the  ther¬ 
mal  shock  (see  Figure  8).  The  damaged  areas  have  no  plated  tungsten  traces.  That  means 
that  the  damage  occurred  after  the  current  was  shut  off.  Therefore,  for  future  runs,  provisions 
will  be  made  not  only  to  increase  the  current  slowly,  but  also  to  decrease  the  current  very 
gradually. 

During  the  first  runs  water  vapor  and  oxygen  content  were  monitored.  The  water  vapor  was 
continuously  sampled  by  a  Beckman  Electrolytic  Hygrometer.  The  initial  water  vapor  content 
in  the  loop  was  well  above  1000  ppm.  It  was  lowered  first  by  purging,  then  by  circulating 
through  the  helium  purification  system.  The  best  value  which  could  be  obtained  was  48  ppm. 

The  water  vapor  of  the  helium,  which  comes  from  the  bottles,  is  11  ppm.  The  additional 
water  vapor  content  in  the  system  probably  comes  from  the  stainless  steel  walls  in  the  system. 
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Figure  1.  Final  form  of  the  schematic  of  the  Allison  closed-loop  MPD  diagnostic  device 
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1.  Charcoal  filter 

2.  NaK  bubbler 

3.  Magnet 

4.  Cesium  tank 

5.  Vacuum  pump 


Figure  2.  A  portion 
of  the  system  located 
inside  the  building— 
during  assembly 


1.  Helium  compressor  3.  Cesium  tank 

2.  Filter  4.  Helium  bottles 


Figure  3.  The  part  of  the  system  located  outside  the  building 
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He  -  Cs  separator 


3.  Second  aftercooler 

4.  Third  aftercooler 


1.  Primary  condenser 

2.  First  aftercooler 


Figure  4.  Cooling  system 


Figure  5.  Control  panel  for  the  system 
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1.  Heater 


2.  Magnet 


Observation  hole 


Figure  6.  Heater  and  magnet  assembly 

The  oxygen  content  was  continuously  sampled  by  an  oxygen  analyzer  (Loekwood  and  McLorie). 
The  best  value  obtained  was  1.  75  ppm.  The  content  of  the  bottle  helium  was  1.  5  ppm. 


The  oxygen  content  should  not  eause  any  trouble;  however,  the  water  vapor  content  should  be 
decreased.  Efforts  are  being  made  to  do  this. 


Figure  9  shows  the  MPD  device  with  the  stainless  steel  insert  during  operation  at  1350°K. 


A  summary  of  the  first  runs  follows: 


1.  The  minimum  mass  flow  for  safe  heater  operation  has  been  determined  to  be  3  gr/see. 

2.  The  heater  operated  satisfactorily  to  at  least  1400°K  gas  temperature. 

3.  The  tantalum  sample  piece  did  not  show  any  wear  at  the  described  operation  conditions. 
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Figure  7.  Tungsten  coil  of  high  temperature 
resistance  heater  before  the  first  run 


Figure  8.  Tungsten  coil  of  high  temperature 
resistance  heater  after  one  hour  running  time 


Discharge  tube 

for  impurity  determination 


Magnetic  coil 


Stainless  steel 
insert 


Figure  9.  MPD-device  with  stainless  steel  insert  during  operation  at  1350°K 
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IV.  DIAGNOSTIC  METHODS 


Spectroscopic  Technique  for  Temperature  or  Pressure 
Measurement  in  Plasmas 


INTRODUCTION 

Spectroscopic  methods  for  measurement  of  temperature  or  pressure  have  the  advantage  that  they 
do  not  change  the  conditions  of  the  plasma.  They  are,  therefore,  very  useful  for  application 
on  small-volume  plasma  devices  (e.g. ,  MHD  power  generators  and  thermionic  energy  conver¬ 
ters). 


Most  of  the  spectroscopic  temperature  measurement  methods  apply,  in  one  way  or  the  other, 
the  basic  relation  for  the  intensity  of  spectral  lines.  Assuming  optically  thin  plasma  and  Boltz¬ 
mann  distribution  of  the  excited  states  of  the  atom,  the  intensity  of  a  spectral  line  resulting 
from  a  transition  m-»n  is  given  by 
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m 

n 


h  V 


m 

n 


n  (P,  T) 
u  (T) 


-E^/kT 
,  m  £ 


where 


(1) 


A  ^  =  transition  probability 

gj„  =  statistical  weight 

h  =  Planck's  constant 

V  =  line  frequency 
n 

n  =  n(P,  T)  =  atom  density 
u  =  u(T)  =  partition  function 
Ej^  =  excitation  energy 
k  =  Boltzmann  constant 
T  =  Temperature  in  “K 
£  =  geometric  length 


By  means  of  Equation  (1),  the  temperature  can  be  determined  from  the  experimentally  observed 
absolute  line  intensity  when  the  transition  probability  and  the  atom  density  are  known.  A  short¬ 
coming  of  this  method  is  that  the  photographic  emulsion  has  to  be  calibrated. 
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The  requirement  for  an  absolute  calibration  can  be  avoided  by  determining  the  temperature 

1 4^ 

from  the  intensity  ratio  of  different  spectral  lines  (Ornstein  ).  For  the  transitions  m-«n  and 
p-»q,  this  intensity  ratio  is  given  by 
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-  (E^-EJ/kT 


(2) 


In  this  case,  only  a  relative  calibration  is  necessary.  Furthermore,  only  relative  transition 
probabilities  have  to  be  known,  for  which,  in  general,  more  accurate  values  are  available  than 
for  absolute  transition  probabilities. 


Relative  temperature  measurement  is  possible,  finally,  without  knowledge  of  transition  proba- 
bilities  (Hoermann  and  Maecker  ).  In  this  method,  use  is  made  of  the  relation  for  the  intensity 
ratio  of  one  and  the  same  spectral  line  at  different  temperatures: 


(T2) 


n  (P,  Tj)  U  (Tj)  - 
n  (P,  T2)  U  (Tj) 


(3) 


Equation  (3)  allows  measurement  of  relative  spatial  temperature  changes  and  absolute  measure¬ 
ment  of  temperature  profiles,  assuming  the  temperature  is  known  at  one  point.  It  is  also  recog¬ 
nized  that  the  sensitivity  of  the  emulsion  need  not  be  known  because  intensities  of  the  same 
wavelength  are  compared. 


If  in  the  plasma  under  consideration  a  spectral  line  has  maximum  intensity  in  a  region  that  is 
known  not  to  contain  the  maximum  temperature,  this  line  intensity  maximum  is  due  to  the  fact 
that,  in  this  case,  the  maximum  radiation  intensity  is  already  reached  at  a  temperature  below 
the  temperature  maximum.  The  temperature  corresponding  to  this  maximum  of  radiation  in¬ 
tensity  can  be  calculated  by  the  method  of  Fowler  and  Milne*  and  used  as  a  reference  tempera- 

5  6 

ture  in  Equation  (3).  Larenz  and  Bartels  used  this  method  to  determine  absolute  temperature 
profiles  in  arcs. 

In  the  following  paragraphs,  a  similar  method  is  described  which,  however,  does  not  require 
the  detection  of  the  maximum  of  intensity  versus  temperature  curve  of  a  spectrum  line.  This 
method  can  be  used  for  plasma  temperature  determination  when  the  total  pressure  of  the  plasma 
is  known,  or  for  pressure  determination  when  the  temperature  is  known. 


♦Superscripts  denote  references  presented  in  Section  VI. 


14 


METHOD 


The  derivation  of  the  basic  equations  for  the  new  method  is  carried  through  for  a  plasma  con¬ 
sisting  of  two  different  gas  species.  As  an  application  in  particular,  an  alkali-seeded  noble 
gas  plasma  is  considered.  For  the  calculation  of  the  particle  densities  in  the  plasma,  the 
validity  of  the  equation  of  state  for  perfect  gases  and  of  Saha's  equation  is  assumed.  In  the 
latter,  a  Boltzmann  distribution  of  the  excited  atomic  states  also  is  assumed.  It  is  remarkable, 
however,  that  these  assumptions  place  no  restriction  on  the  method  itself  which  can  also  be 
extended  to  plasmas  deviating  strongly  from  thermodynamic  equilibrium. 

Let  the  densities  of  the  neutrals  and  ions  be  designated  by  n^,  n^^,^  for  the  first  gas  species, 
and  by  ng,  ng.^  for  the  second  gas  species.  The  total  electron  density  in  the  plasma  is  n.. 

With  these  designations  the  Saha  equations  describing  the  ionization  equilibrium  in  the  two- 
gas  plasma  are: 


"a+"- 


'  4>k>  4>a  =  2 


/2jr^\3/2 

V 


where 


=  effective  ionization  energy  of  the  neutrals  A 
Wg  =  effective  ionization  energy  of  the  neutrals  B 


^A  =  E  gAs  e 

s  =  0 


=  E  gAs  ® 


/xAs\ 

\kT) 


neutrals  A 


partition  functions  of  the 


ions  A+ 
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neutrals  B 

s  =  0 

/x 
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- 

partition  function  of  the 

CO 

.  -ft 

ub+  "  £ 

g^  e  V  ^ 

®Bs 

T/ 

ions  B+ 

s  =  0 

^As’ 


®Bs' 


ggg  =  statistical  weights  of  the  intrinsic  states  of  the 
neutrals  and  ions 


^As*  *Bs* 


"Bs 


=  excitation  energies  of  the  intrinsic 
states  of  the  neutrals  and  ions 


The  considered  ionization  equilibrium  refers  to  a  temperature  range  and  mixing  proportions 
for  which  double  and  higher  multi -ionizations  can  be  neglected  quantitatively.  Equations  (4-5) 
have  to  be  supplemented  by  the  equation  of  state: 

P  =  (n^  t  n^^  +  ng  +  ng^  +  n  ,)kT  (equal  temperature  of  the  components  assumed)  (7) 

the  condition  for  electrical  neutrality: 

=  "a+  ^  "B+  <8) 

and  the  definition  equation  for  the  mixing  proportion: 

Ha  +  Ha  + 

a  =  _2 - ill  (9) 

+  ‘^B  + 


From  Equations  (4)  through  (9),  the  particle  densities  can  be  calculated  as  function  of  tempera¬ 
ture,  T,  total  pressure,  P,  and  mixing  proportion,  a  .  It  results  that  the  electron  density  is 
given  as  the  positive  real  root  of  the  cubic 


3^2 
n  +  n  . 


•^A  +  ^B  +  -  ^ 


1+  a 


1+  a  ■  kT 


(10) 


16 


The  other  particle  densities  can  be  reduced  to  the  electron  density: 


The  further  relations  are  equal  in  form  for  both  neutral  species.  For  convenience,  the  in¬ 
dices  A  and  B  will,  therefore,  be  omitted. 


By  means  of  Equations  (10),  (13),  and  (14),  the  following  expression,  consisting  of  particle 
densities  and  partition  functions  of  one  of  the  neutral  species,  can  be  calculated  as  function  of 
pressure  P  and  temperature  Tj,  with  a  as  parameter; 


b  =  log 


■n(P,  Ti)  u(T2)' 

n(P,  T2)  u(Ti)J' 


T2  = 


1-aTi 


Physically,  the  expression  b  refers  to  two  points  in  the  plasma  with  different  temperatures, 
Tj  and  T2.  The  pressure,  P,  is  presumed  constant.  This  is  valid  for  many  applications 
where  exterior  force  fields  are  negligible  [otherwise,  the  substitutions  (P,  Tj) -»(Pi ,  Tj ), 
(P,  T2)  —  (P2.  T2)  have  to  be  made  in  Equation  (15)]  .  The  expression  b  and  the  temperature 
parameter, 

1  1 


are  experimental  observables  and  can  be  measured  by  spectroscopic  methods.  Thus,  by  in¬ 
serting  the  experimentally  determined  a-  and  b-  values  into  Equation  (15),  an  implicit  equation 
in  P  and  T^  is  obtained.  From  this,  either  Ti,  when  P  is  known,  or  P,  when  T^  is  known, 
can  be  calculated  numerically. 


Let  us  now  explain  how  the  fundamental  expression  b  and  the  temperature  parameter  a  are  de¬ 
termined  spectroscopically. 
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Assuming  a  distribution  according  to  Boltzmann  statistics,  the  number  density  of  atoms  in  ex¬ 
cited  state  m  is  given  by 


n 


m 


n(P,  T) 
u  (T) 


-Em/kT 

gm  e 


(17) 


For  the  ratio  of  number  densities  of  excited  atoms  in  the  state  m  in  the  same  plasma  at  dif¬ 
ferent  temperatures,  Tj  and  T2,  it  follows  that 


^m  1 

nm(Ti)  .n(P,  Tj)  u(T2)  '  ^  ' 

Ui 

T2) 

(18) 

nm  ’  n(P.  T2)  '  ufTj) 


According  to  Equations  (D  and  (18).  the  intensity  ratio  of  two  spectral  lines  measured  at  two 
points  with  different  temperatures,  Tj  and  T2.  in  the  same  plasma  and  resulting  from  one  and 
the  same  transition.  m-*n,  is  given  by 

1"(T,)  n(P.  T,|  ,  (T,|  ■%.  (19 

(T2)  ’  n  (P,  T2)  u  (Ti)  ® 


Dividing  Equation  (19)  by  the  corresponding  equation  for  the  transition  p-*q,  the  following  rela¬ 
tion  is  obtained: 


(20) 


By  introducing  the  temperature  parameter  a  [see  Equation  (16)3and  the  decadic  logarithms  of 
the  intensity  ratios: 


AY 


m 


log 


C  <Ti) 

IJT  <T2) 


log 


(Ti) 

(T2) 


(21) 


Equation  (20)  can  be  rewritten  in  the  form: 


•  a 


k  In  10 


AY*^  -  AYP 

n  q 


(22) 
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This  relation  permits  the  temperature  parameter  a  to  be  determined  spectroscopically  by 
measuring  the  decadic  logarithms  of  the  intensity  ratios  defined  in  Equation  (21)  for  two  spectrum 
lines.  Similarly,  by  introducing  the  fundamental  expression  b  [see  Equation  (15)]  ,  Equation 
(19)  can  be  rewritten  in  the  form: 


b  = 


k  In  10 


(23) 


This  relation,  Equation  (23),  shows,  finally,  that  the  fundamental  expression  b  can  be  deter¬ 
mined  experimentally  from  the  measured  -  and  a-  value.  The  demonstrated  method 

works  without  knowing  transition  probabilities.  According  to  Equations  (22)  and  (23),  only  the 
two  excitation  energies,  E^  and  Ep,  have  to  be  known.  The  compares  intensities  at 

the  same  wavelength;  therefore,  the  dependence  of  sensitivity  from  the  wavelength  does  not 

have  to  be  known. 


Usually,  in  most  plasmas  the  change  in  photographic  density  is  much  greater  due  to  the  change 
in  sensitivity  of  the  photographic  plate  rather  than  to  the  change  in  radiation  from  the  plasma- 
using,  for  example,  the  line  ratio  technique  (Equation  2).  It  is,  therefore,  very  desirable  to 

remove  the  source  of  error. 


As  already  mentioned,  by  inserting  the  experimentally  obtained  a-  and  b-  value  into  Equation 
(15),  either  Tj  can  be  calculated  if  P  is  known,  or  P  can  be  calculated  if  T i  is  known.  In 
practical  applications,  however,  it  is  suitable  to  make  use  of  a  graphic  representation.  For 
this,  the  theoretical  curves  representing  the  fundamental  expression 


are  drawn  as  function  of  Ti  and  P,  respectively  (abscissa),  where,  in  the  first  case,  a  and  P  and 
in  the  second  case,  a  and  Tj  are  treated  as  parameters.  From  these  two  parametric  curve 
representations,  the  unknown  temperature,  Tj,  or  pressure,  P,  respectively,  are  read  as  fol¬ 
lows.  The  experimental  a-  value  and  the  second  parameter,  P  or  Tj,  respectively,  designate 
the  theoretical  b-  curve  that  applies.  The  experimental  b-  value  gives  a  definite  point  on  this 
curve.  The  abscissa  of  this  point  represents  the  unknown  observable. 
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As  an  application,  the  theoretical  b-  value  has  been  calculated  as  a  function  of  temperature 
Tj  =  1000  to  6000°K  for  a  cesium-seeded  helium  plasma  —  with  a  pressure  of  P  =  1  atmosphere, 
the  seeding  ratio  a  =  10"^,  and  a  number  of  a-  values  as  parameters.  The  ionization  potentials 
of  cesium  and  helium  were  taken  as 

=  3,893  eV,W„  =  24.580  eV  (25) 

C  s  1*6 

The  depression  of  the  ionization  potentials  due  to  the  electrostatic  microfield  and  the  Debye 

7 

polarization  representing  under  the  given  conditions  a  small  correction 

AW  =  0.67  X  10‘®  nJ +  0.37  X  lO"”^  eV,  (26) 

was  considered  by  means  of  a  successive  approximation. 

The  numerical  results  are  represented  graphically  in  Figure  10  for  a  helium-cesium  mixture 
with  a  seeding  ratio  a  =  10'^.  Curves  for  pure  cesium  and  pure  helium  and  other  seeding  ratios 
will  be  included  in  a  later  report. 


ERROR  CONSIDERATION 


The  accuracy  of  the  method  depends  noticeably  on  the  compared  temperatures  Tj,  T2  in  the 

1  1 

plasma  — i.e.  ,  also  on  the  temperature  parameter  a  =  — ^  -  —  — and  the  pressure  P,  With  re¬ 
spect  to  temperature  measurement,  this  is  recognized  easily  from  the  theoretical  b-  curves 
(Figure  10)  by  estimating  graphically  the  error  ATj  due  to  certain  assumed  experimental 
errors,  Aa  and  Ab,  over  the  whole  range  of  Tj,  a,  and  b.  To  see  how  the  error  in  the  meas¬ 
uring  result  depends  (in  principle)  on  the  respective  values  of  Tj,  P,  and  a,  an  analytical 
representation  is  given  for  a  simple  plasma  consisting  of  one  kind  of  atomic  and  ionic  species. 
In  this  case,  the  neutral  density  is  given  by  [see  Equations  (4),  (7),  and  (8)]  : 


(27) 
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Electron— Temperature—  K 


Figure  10.  Theoretical  b-curves  for  atomic  cesium  lines  for  a  helium-cesium 
mixture.  Seeding  ratio  a  =  10'^,  Total  pressure  =  1  atmosphere 
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and  the  fundamental  expression  b  by  [see  Equation  (15)]  : 


The  maximum  errors  in  a  and  b  due  to  the  measuring  errors  of  AYP  and  E  E  are 

n  q  m’  p 

given  by  the  relation  [see  Equations  (22)  and  (23)]  : 


Evaluating  the  total  differential  of  Equation  (28)  with  respect  to  a,  b,  Tj,  and  replacing  the 
differentials  by  the  corresponding  erros,  the  theoretical  relation  between  the  errors  Aa,  Ab 
andATj  and  AP  is  obtained.  From  this  relation,  the  maximum  relative  error  ATj/Tj  in 
temperature  and  the  maximum  relative  error  AP/P  in  pressure  are  derived  as: 


(31) 


(32) 


22 


(33) 
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Equations  (31)  and  (32)  give  the  relative  error  ATj /Tj  of  the  temperature  measurement  and  the 
relative  error  AP/P  of  the  pressure  measurement  as  function  of  P,  Tj.  a  and  the  corresponding 
measuring  errors  {  Aa,  Ab,  AP  }  and  |  Aa.  Ab.  ATj  },  respectively.  Thus,  by  estimating  the 
measuring  errors,  the  error  in  the  measuring  result  (i.e.,  the  exactness  of  the  method)  in  the 
respective  region  of  P,  T,  and  a  can  be  precalculated. 

Numerical  investigation  of  the  error  relations.  Equations  (31)  and  (32),  show  that  the  method 
can  be  used  with  a  considerable  degree  of  accuracy  [(ATj/Tj)  100  <10%;  (AP/P)  100  <  10%]  ,  if 
the  plasma  temperature  or  plasma  pressure  lies  in  that  temperature  or  pressure  region,  re¬ 
spectively,  where  the  theoretical  b  versus  Tj  curve  (temperature  measurement)  or  theoretical 
b  versus  P  curve  (pressure  measurement)  is  changing  rapidly.  Furthermore,  it  is  found  that 
the  method  is  the  more  accurate  the  larger  the  experimentally  observed  a-  parameter  is  — i.e,  , 
the  more  different  the  temperatures  of  the  compared  spectral  lines  are.  However,  it  has  to  be 
considered  that  at  very  high  a-  values,  the  plasma  pressure  can  only  be  assumed  to  be  approxi¬ 
mately  constant.  In  this  connection,  it  is  remarkable  that  the  temperature  measurement  is  not 
very  sensitive  to  small  uncertainties  in  the  plasma  pressure. 

With  regard  to  a  1%  cesium  seeded  helium  plasma  at  atmospheric  pressure,  the  method  is  re¬ 
stricted  to  a  temperature  region  between  Ti  S  2000  to  5000»K.  This  is  exactly  the  region  of 
tne  electron  temperature  within  which  MPD  energy  conversion  becomes  of  practical  interest. 
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Calculation  and  Measurement  of  the  Radial  Distribution 
of  Line  Broadening  in  a  Cesium  Plasma 

In  MPD  power  generation  devices  as  well  as  in  thermionic  energy  converters,  the  plasma 
properties  change  (in  general)  rapidly  within  the  electrode  sheaths  and  boundary  layers. 
Therefore,  a  method  for  the  determination  of  electron  density  in  nonuniform  plasmas  is  de¬ 
sirable,  One  such  method  for  determination  of  the  electron  density  is  based  on  the  line  broad¬ 
ening  caused  by  the  Stark  effect.  In  the  case  of  cesium,  which  has  been  used  in  both  .MPD 
power  generators  and  thermionic  converters,  this  method  has  been  applied  successfully  by 
several  investigators,*'*  Theoretical  calculations  of  the  line  width  for  cesium  have  been  made 
by  Griem,**  Agnew,  and  Summers.” 

In  the  interpretation  of  spectrum  it  must  be  recognized  that  the  observed  radiation  comes  from 
different  plasma  layers,  where  different  electron  densities  are  present.  Therefore,  the  line 
width  observed  on  the  photographic  plate  may  be  a  mixture  of  a  large  number  of  differing  line 
widths.  In  the  case  of  a  cylindrical  plasma  configuration  it  is  possible  to  separate  this  mix¬ 
ture.  This  procedure  is  described  later  in  this  subsection. 

In  addition  to  determining  the  true  line  width  rather  than  an  averaged  value,  this  procedure 
allows  several  line  widths  corresponding  to  the  electron  density  profile  of  the  discharge  to  be 
obtained  through  the  analysis  of  a  single  line.  This  is  extremely  important  in  applications 
which  allow  a  short  running  time  and  necessitate  a  maximum  of  information  for  each  single 
exposure, 

MATHEMATICAL  METHOD 

The  problem  of  determining  the  true  intensity  is  due  to  the  fact  that  a  cylindrical  gas  discharge 
usually  has  a  temperature  and  density  profile  and,  therefore,  does  not  radiate  uniformly.  This 
problem  has  been  approached  and  solved  for  the  case  of  continuous  and  line  radiation,  neglect¬ 
ing  broadening,  by  Hermann,*  Maecker,*  and  Pierce,”  The  resulting  solutions  have  been  used 
by  many  other  investigators. 

To  obtain  the  true  distribution  of  intensity  within  the  discharge,  it  is  necessary  to  consider 
the  Abel  integral  equation.  This  equation  is  based  on  the  transformation  of  intensity  from 
cartesian  to  cylindrical  coordinates; 
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dy  r 

In  differential  form,  this  equation  becomes  = 

dr 

The  intensity  at  any  point  x.  y  is  related  to  the  intensity  density  i(r)  by 


d  Ij^^  y  =  i(r)  dy 

r  dr 

It  follows  that:  d  I„  ,,  =  i(r)  /— s - rr 

X.  >  ^^2  . 


Integration  of  the  expression  between  the  limits  r  =  x  and  r  =  R  results  in: 


I  (x)  =  2 


R 

/ 

r=x 


d  I 


X.  y 


A  factor  of  2  has  been  added  in  the  last  expression  to  include  both  halves  of  the  cylinder. 


Figure  11  illustrates  the  usual  approach  for  solution  of  this  equation  by  numerical  integration. 
The  notation  is  after  Pierce.**  In  Figure  1 1,  the  emitting  volume  is  assumed  to  be  a  uniformly 
thick  slab  with  cylindrical  geometry.  The  density  of  the  emitter  is  assumed  to  be  a  step  function 
of  the  radius  rather  than  a  continuous  function.  The  intensity  density,  assuming  a  thickness  of 
unity,  of  the  segment  (j,  k)  becomes; 

Ij 

i  =  - ^ - ,  where  L  is  the  observed  intensity- 

(j,  k)  2AAj_  ^  •’ 

The  notation  indicates  an  axial  location  j  and  radial  location  k.  Since  the  discharge  may  be 
considered  as  axisymmetrical,  the  intensity  of  the  outer  ring  (Figure  11)  becomes 

*k  ■  *(j,  k)  ■  2AAj  y 

The  intensity  density  of  the  second  ring  referenced  from  the  outer  diameter  (Figure  11)  may  be 
found  by  subtraction  of  the  radiated  energy  of  the  outer  ring  from  the  observed  intensity  at 
position  j-1.  Therefore, 

^j-1  "^^k^^j-l,  k 

*k-i  = 

k-1 
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RADIAL  SEGMENTS  OF  CYLINDRICAL  DISCHARGE 


STIGMATIC  LINE 


Figure  11.  Numerical  integration. 


As  progression  is  made  toward  the  center  of  the  discharge,  the  intensity  density  within  each 
incremental  ring  may  be  found.  The  intensity  densities  of  each  radial  increment  multiplied  by 
the  respective  transverse  areas  are  used  to  calculate  the  intensity  density  in  the  next  smaller 
ring.  By  selecting  small  increments  the  procedure  provides  a  profile  which  approaches  the 

true  radial  distribution  of  intensity  density.  The  scries  expression  for  this  computation  be- 
comes: 


^k-n 


Ij.n  -  k-(n-l)  ~  ^^^j-n,  k-(n-2)-- • 


-  k 


k-n 


where  n  =  1,  2,  3,  4, 
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Each  incremental  area  as  illustrated  in  Figure  11  may  be  computed  by 

"  2  k  j,k  ^  ■  ^k-1^  k-1  '^^j,  k-P  + 

^^j.k  ■  ^j.  k-P  ■  ^j-1  <yj-l,  k  "yj-l,  k-l^] 

where 

Xi 

6=  cos"^  — 
r; 


and 


y  =  rj^  sinfl 

The  foregoing  procedure  converts  an  experimentally  obtained  intensity-versus-x-distribution 
in  an  intensity-density -versus-r-distribution.  This  procedure  can  be  applied  either  on  con¬ 
tinuous  radiation  or  on  line  radiation  without  broadening. 

In  order  to  take  into  account  the  broadening,  the  procedure  becomes  more  complicated.  Figure 
12  is  a  schematic  of  the  computational  procedure.  Figure  12a  is  a  sketch  of  a  typical  line  of  an 
observed  stigmatic  spectrum.  The  line  is  normally  more  broadened  in  the  inside  than  in  the 
outer  zones  where  the  electron  density  gradually  goes  to  zero.  This,  however,  is  not  always 
true— e.  g,  ,  in  high  temperature  plasmas  lines  with  a  torus-like  shape  have  been  observed. 

The  line,  regardless  of  its  shape,  is  divided  in  ten  zones  as  indicated  in  Figure  12a.  After 
scanning  with  a  microdensitometer,  a  family  of  traces  of  photographic  density  versus  d  is  ob¬ 
tained  (Figure  12b). 

In  the  next  step  the  photographic  density  is  converted  into  intensity  by  the  relation: 

D  =  spectrum  density 
D  -  y  log  I  I  =  intensity 

y  =  gamma  function  of  photographic  plate 
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Since  this  relation  holds  only  in  the  density  range  0.  8-1.  8,  in  this  application  it  is  necessary 
to  use  Seidel's  transformation. 


to  convert  Figure  12b  to  Figure  12c. 

T  is  related  to  D  by; 

1 

D  =  log  - 

The  family  of  curves  obtained  in  Figure  12c  is  now  an  intensity-versus-X-distribution  with  x 
as  parameter. 

The  next  step  is  to  transform  this  into  a  family  of  curves  with  x  as  ordinate  and  X  as  parameter. 
This  change  of  variables  is  accomplished  by  reading  in  all  the  x's  for  a  given  X.  This  procedure 
is  a  common  computer  problem  and,  therefore,  will  not  be  described  further  here.  Now,  the 
conventional  Abel  transformation,  already  described,  is  applied  in  going  from  Figure  12d  to 
Figure  12e.  A  family  of  curves  is  now  obtained  which  gives  an  intensity-density-versus-r 
distribution  with  X  as  parameter. 

To  go  from  Figure  12e  to  Figure  12f,  again  a  change  of  variables  has  to  be  made.  For  a  cer¬ 
tain  r  all  the  different  X  values  are  read  in  and  a  family  of  curves,  intensity-density-versus-X. 
is  obtained.  Thus,  the  mixture  of  all  the  different  line  shapes,  which  originally  existed,  is 
separated  into  its  components.  The  last  step  is  to  plot  the  half  width  of  the  lines  versus  r 
(Figure  12g).  Using  the  results  of  Griem’®or  Agnew  and  Summers,"  this  plot  can  be  converted 
into  an  electron  density  profile. 

APPLICATION 

The  numerical  procedure  has  been  applied  to  line  broadening  measurements  of  a  cesium  diode, 
the  spectrum  of  which  has  been  taken  with  an  Ebert  3.4-meter  grating  spectrograph.  Data  have 
also  been  obtained  with  respect  to  line  broadening  of  helium  in  an  arc  discharge.  The  true  dis¬ 
tribution  of  half  width  (width  at  one-half  maximum  intensity)  is  compared  with  the  distribution  ■ 
obtained  from  the  observed  stigmatic  spectrum  (Figure  13). 
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Figure  13.  Radial  distribution  of  half-width  k  =  6629  A. 

In  order  to  ensure  independence  of  nonequilibrium  effects,  the  fundamental  series  line,  GG29  A, 
has  been  isolated  for  the  analysis.  Note  from  Figure  13  that  the  observed  distribution  tends  to 
give  an  average  of  the  true  distribution. 

The  distribution  of  half  width  has  been  converted  to  the  radial  distribution  of  electron  density 
according  to  the  relation  presented  by  Agnew  and  Summers, 

Xg  =  2’^'^  X  1, 5  Ng  X  lO'^^ 

For  the  fundamental  series  transitions.  5D-nF  for  5  ^n  ^8,  the  line  broadening  has  been  esti¬ 
mated  to  be  independent  of  electron  temperature  over  a  temperature  range  of  2000°K  to  10,  000°K. 
Since  n  =  8  for  the  6629  A  line,  the  electron  density  is  expressed  as; 

1  c  electrons 
Np  =  -  X  10^5  - 5 - 
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By  applying  this  equation,  the  radial  distribution  of  electron  density  as  shown  in  Figure  14  is 
obtained.  The  electron  density  is  low  in  the  center  of  the  discharge  where  the  temperature  is 
maximum.  For  this  type  of  discharge,  the  partial  pressure  (P^  =  NgkTg)  tends  to  remain 
constant.  As  the  temperature  decreases,  tiie  electron  density  increases  until  a  point  is 
reached  where  the  temperature  is  siiffieiently  low  that  ionization  becomes  negligible.  Thus, 
the  eleetron  density  starts  to  decrease,  Tlie  rise  of  electron  density  at  the  outer  radial  zones 
is  probably  duo  to  the  back  currents  in  the  walls  of  the  plasma  container. 

This  method  allows  a  detailed  investigation  of  radial  distribution  of  electron  density.  Because 
the  stigmatic  line  spectrum  docs  not  allow  precise  determination  of  the  location  of  the  electron 
distribution  within  the  discharge,  the  method  is  used  in  conjunction  with  the  Spectrum  Picture 
Method  discussed  in  the  hhrst  Quarterly  Technical  Summary  Report,  EDR  3155. 
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V.  THEORETICAL  PLASMA  PROPERTIES 


The  computations  of  the  properties  of  cesium-seeded  helium  plasma.  Tables  I  and  II  and 
Figure  15,  were  carried  through  by  consideration  of  the  ionization  of  cesium  and  helium.  The 
theoretical  relations  used  are  given  in  Equations  (1),  (6),  (10)  through  (14),  and  (26).  The 
ionization  potentials  were  corrected  for  the  depression  by  the  electrostatic  microfield  and  the 
Debye  polarization.  The  main  assumptions  are  distribution  of  the  excited  states  according  to 
Boltzmann  statistics  and  existence  of  local  thermodynamic  equilibrium.  The  variables  con¬ 
sidered  were: 


Temperature,  T 
Seeding  Ratio,  o 
Pressure,  P 


1000  to  15,  OOO-’K 
10"^ 

1  atmosphere 
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Particle  densities. 
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Table  II. 

Atomic  partition  functions. 
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